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ABSTRACT 

We analyze multi-band photometry of a subsample of low mass stars in the 
associations Ori OBla and lb discovered during the CIDA Orion Variability 
Survey, which have ages of 7 - 10 Myr and 3-5 Myr, respectively. We obtained 
UBVRcIc photometry at Mt Hopkins for 6 Classical T Tauri stars (CTTS) and 
26 Weak T Tauri stars (WTTS) in Ori OBla, and for 21 CTTS and 2 WTTS in 
Ori OBlb. We also obtained L band photometry for 14 CTTS at Mt. Hopkins, 
and 10/im and 18/im photometry with OSCIR at Gemini for 6 CTTS; of these, all 
6 were detected at lO/im while only one was detected at 18/im. We estimate mass 
accretion rates from the excess luminosity at U, and find that they are consistent 
with determinations for a number of other associations, with or without high 
mass star formation. The observed decrease of mass accretion rate with age is 
quahtatively consistent with predictions of viscous evolution of accretion disks, 
although other factors can also play a role in slowing accretion rates. We compare 
the excesses over photospheric fluxes in H — K, and K — L, and K — N with 
the younger sample of Taurus and find an overall decrease of disk emission from 
Taurus to Ori OBlb to Ori OBla. This decrease implies that significant grain 
growth and settling towards the midplane has taken place in the inner disks of 
Ori OBI. We compare the SED of the star detected at both lO/im and IS/xm 
with disk models for similar stellar and accretion parameters. We find that the 
low < 18/im fluxes of this Ori OBlb star cannot be due to the smaller disk 
radius expected from viscous evolution in the presence of the FUV radiation 
fields from the OB stars in the association. Instead, we find that the disk of this 
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star is essentially a flat disk, with little if any flaring, indicating a a significant 
degree of dust settling towards the midplane, as expected from dust evolution in 
protoplanetary disks. 

Subject headings: accretion — stars: pre-main sequence — stars: formation — 
infrared: stars — techniques: photometric 

1. Introduction 

Most star formation studies are based on a few nearby associations without massive 
stars, as Taurus and Chameleon or in highly populated clusters with significant high mass 
formation, like the Orion Nebula Cluster (ONC). In recent years, several teams have begun 
systematic studies of the populations of the OB associations. These associations cover large 
volumes in the sky and include populations covering the entire mass range. OB associations 
have been subject of studies with Hipparcos (Brown et al. 1997; de Zeeuw et al. 1999), 
which has helped pinpoint their distances and mean ages. In a number of associations, the 
populations display a gradient of ages, from ~ 1 Myr to a few ~ 10 Myr, which has been 
interpreted as the result of triggered star formation (Blaauw 1964). This range of ages makes 
the OB associations the most suitable laboratory for studies of protoplanetary disk evolution, 
because it encompasses the time span in which giant planets arc expected to form (Pollack 
et al. 1996; Alibert et al. 2004), and significant solid evolution is expected to take place, as 
indicated by meteoritic evidence (Podoseck & Cassen 1996; Wood 2004). 

We are carrying out a photometric and spectroscopic survey of ~ 128 degrees in the 
Orion OBI association in order to identify the low mass population. Initial results of this 
survey have been presented by Briceno et al. (2001, 2004 = Paper I). In Paper I, we 
identified 197 new low mass members in the subassociations Ori OBla and lb, spanning the 
mass range from ~ 0.2 to 1.4 Mq. We found ages of ~ 3 - 5 Myr for Ori OBlb and ~ 7 - 
10 Myr for Ori OBla, confirming age determinations based on the OB stars (Blaauw 1964; 
Warren & Hesser 1977, 1978; Brown et al. 1994). We also found that Ori OBlb could be 
identified with a population associated with a ring of low density gas and dust, probably 
formed by a supernova event, with most of the Classical T Tauri stars (i.e., stars accreting 
from circumstcUar disks; CTTS) near the ring, similar to the case of A Ori (Dolan & Mathieu 
1999, 2001, 2002). In contrast, Ori OBla is in a region devoid of gas and dust, in agreement 
with its older age. We also found that the number of CTTS decreases sharply between Ori 
OBlb and la, indicating the few disks remain by ~ 10 Myr, in agreement with findings 
from young clusters, some containing high mass stars, and low mass associations as TW 
Hya (Haisch, Lada, & Lada 2003; Hillenbrand, Carpenter, and Meyer 2004, in preparation; 
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Sicilia-Aguilar et al. 2004, in preparation; MuzeroUe et al. 2000, 2001). 

In this work, we present multiwavelength photometry of a subset of the stars in Paper 
I, including stars in Ori OB la, lb and Ic. In particular, we have obtained U photometry 
of a significant number of CTTS in these subassociations, from which we have determined 
mass accretion rates for comparison with disk evolution theories. We have also obtained L 
band and mid-infrared (IR) measurements at 10/xm and 18/im of a small subset of stars in 
Ori OBI, which has allowed us to make assessments of the state of the dust in their disks 
in comparison with younger populations. In §2 we present the photometry, and in §3, we 
examine the properties of the disks in Ori OBla and lb in comparison with the Taurus 
populations. We discuss the implications of these results for disk evolution in §4. 

2. Observations and Data Analysis 

The targets were selected from the sample of the Orion variability survey presented 
in Paper 1, which were originally found on the basis of variability in the V band using 
the QuEST camera on the Im Schmidt telescope on the Venezuela National Observatory. 
Follow-up spectroscopy using the FAST Spectrograph (Fabricant et al. 1998) on the 1.5m 
telescope of the Fred Lawrence Whipple Observatory (FLWO) at Mt. Hopkins, and the 
Hydra multi-fiber spectrograph (Harden & Armandroff 1995) on the WYIN/3.5m telescope 
at Kitt Peak, confirmed their membership to Ori OHl (cf. Paper I). 

2.1. Optical photometry 

We used the 4SH00TER CCD array on the 1.2m telescope of the Fred L. Whipple 
Observatory on Mt. Hopkins during two observing runs to obtain UBVRcIc photometry of 
our hst of stars. The first run was during the nights of November 29 through December 4, 
2002, although observations could be done only in the second half of the night of December 
3 and all night on December 4. The second run was during October 15 - 20, 2003. All six 
nights were clear. The seeing was between 1.2" and 2.5" throughout the observations. The 
4SH00TER camera contains four 2048 x 2048 Loral CCDs separated by 45" and arranged in 
a 2 X 2 grid. After binning 2x2 during readout, the plate scale was 0.67" pixel"^. In order to 
achieve more uniform and consistent measurements, we placed all the target standard stars 
on the same CCD detector (chip 3). 

We obtained BVRfJc photometry for 65 stars in Ori OBla, lb, and Ic, and [/-band 
photometry for 59 of these stars. In particular, out of 36 stars in Ori OBla, we obtained 
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U photometry for all 6 CTTS and 26 WTTS; for the 27 stars in Ori OBlb, we obtained 
U photometry for the 2 WTTS and 21 CTTS; finally, we obtained U photometry for the 2 
CTTS observed in Ori OBlc. 

The basic processing of the FLWO data was done with IRAF^ routines in the standard 
way. For both runs, the [/-band data were flat-fielded using sky fiats taken at dusk/dawn. 
Instrumental magnitudes were obtained using the APPHOT package in IRAF. Though Orion 
is at a moderately low galactic latitude b ~ —20°, our fields are not crowded in the relevant 
magnitude range, so aperture photometry is adequate. We used an aperture radius of 15 
pixels for the source, ~ 4x the typical 3.5 pixel FWHM of our observations, large enough 
to include all the light from the standard stars as well as some observations made under 
mediocre seeing. We used a sky annulus with of 15 and 20 pixels for the inner and outer 
radii, respectively. The instrumental magnitudes were calibrated in the Johnson UBV and 
Cousins RI system with observations of Landolt (1992) standard fields SA92, SA95, SA113 
and PG0231+051 at various air masses every night; each field contains stars with a range 
of colors similar to our target objects. Our overall photometric errors are dominated by the 
transformation uncertainty, which are 0.05 mag in U and 0.02 magnitudes for BV{RI)c- 
The nightly variation in the zero points was 0.15, 0.07, 0.06, 0.04, 0.05 magnitudes for 
UBV{RI)c respectively, suggesting reasonably good photometric stability throughout both 
runs. 

The resultant photometry is listed in Table 1, where stars are identified by their running 
numbers in Paper 1. Two exposures in some bands were obtained for star CVSO 35 in Ori 
OBla and for stars CVSO 146 and 190 in Ori OBlb, as indicated in Table 1. 

2.2. Near-infrared photometry 

We obtained L-band photometry for a subsample of 14 stars using the dual channel 
infrared camera STELlRCam (Tollestrup & Willner 1998) on the 1.2m telescope at FLWO 
during the nights of December 18-19, 2002. STELlRCam consists of two 256 x 256 pixel 
InSb detector arrays. Each is fed from a dichroic mirror that separates wavelengths longer 
and shorter than 1.9//m into two independent imaging channels (red and blue channels) 
for simultaneous observations on the sky. Three separate magnifications can be selected 
by rotating cold lens assemblies into the beam. The magnification is the same in both 
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wavelength channels. For our observations we used the medium lens yielding a scale of 
0.6pixel~^ and a 2/5 x 2/5 field of view, together with the L-Barr filter (3.5/im) in the red 
channel. To convert the instrumental magnitudes of STELlRCam to the standard system, 
stars BD+0.1694, G77-31, Gl 105.5, Gl 406, HD 1160, HD 18881, HD 225023 and HD 40335 
from the list of Elias et al. (1982) were observed each night at various air masses. Both the 
Orion and standard stars were observed with a 3 x 3 square pattern with 20" dither between 
positions. The integration time at each dither position was 10 s (0.1 x 100 coadditions) . 

All data were first linearized using frames obtained at various exposure times and rou- 
tines developed in the Interactive Data Language (IDL) by T. Mcgcath at CfA. We then 
used standard IRAF routines to proceed with the data reduction. Average dark frames were 
constructed from darks taken at the beginning and end of each night's observations. These 
average dark frames were then subtracted from each image. Sky frames were individually 
made for each observation by median-combining the nine unregistered frames in one dither 
pattern. The sky frames were subtracted from each image within the dither set. The high 
and low airmass frames were grouped into pairs, and the low airmass frame in each pair was 
subtracted from the high airmass frame. These subtracted pairs were then normalized and 
the pairs combined using a median statistic to create the flat fleld. 

Typically only the target star was visible in the L-band images, with a FWHM ranging 
from ~ 2.3 to 3 pixels. Thus, wc measured instrumental magnitudes with the IRAF AP- 
PHOT package, using an aperture radius of 6 pixels, and a sky annulus with inner and outer 
radii of 10 and 20 pixels respectively, for both Orion and standard stars. After obtaining 
the zero points and color terms (which were small) we arrived at an L-band la error = 0.07 
mags, mostly dominated by measurement errors. 

The L-band photometry is presented in Table 1, where we also include 2MASS JHKg 
magnitudes for the stars from Paper I. 

2.3. Mid-infrared photometry 

During the nights of December 3, 5 and 9, 2003 we obtained 10//m and 18/im photometry 
of a subsample of 6 stars in Ori OB la and B. This subsample was selected among those CTTS 
with Kg < 11 and with large excesses in the JHKg diagram. The observations were obtained 
in service mode using OSCIR on Gemini North at Mauna Kea. OSCIR is a mid-infrared 
(8 — 25/xm) imager and low/medium- resolution (R=100-1000) spectrograph, which uses a 
128 X 128 pixel Si: As IBC detector optimized for the wavelength range 8 — 25/xm. On Gemini, 
OSCIR has a plate scale of about 0.084 arcsec/pixel and a total field of view in its imaging 
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configuration of 11" x 11". We used the N-wide filter (Ag// = 10.8//m and AA = 5.23//m) 
and the IHW18 filter (Ae// = 18.2/im and AA = 1.65/im). In the rest of this article we will 
refer to these filters as the 10/im and 18/im filters. Typical exposure times for the Orion 
targets were 60s at lO/im and 300s at ISfim. Sky conditions were mostly clear through the 
observations. All the observations were performed using the standard technique of chopping 
and nodding, with a chop throw of 15" in declination. For calibration of the photometry the 
standard star a CMa (Sirius) was observed on all nights at both wavelengths. Flux density 
estimates for Sirius were calculated using the SEDs pubhshed by Cohen et al. (1999). For 
Sirius, integration times were of the order of 6s in both lO/im and IS/im. The data were 
processed through the standard pipeline for OSCIR data. Chop pairs are subtracted: beam 
A - beam B or beam B - beam A for each set, depending on the nod position. Then these 
differences are averaged to form a final image. Because we were dealing with point sources, 
per recommendation of the Gemini staff no fiat fielding was performed. The OSCIR array is 
very fiat over the small number of pixels within the point spread function. Experience shared 
between Gemini and University of Florida is that the errors associated with the photometric 
calibration (~ 10%) are much larger than any fiat field correction. 

The lOfim and 18/xm images showed only one object per frame, when detected at all. 
We used APPHOT to obtain aperture photometry of all targets. Images at lO/im showed 
typical FWHM~ 6 pixels (or ~ 0.5 arcsec), so we used an aperture radius of 15 pixels and 
inner and outer sky annulus radii of 15 and 25 pixels respectively. Because most of the 
Orion stars were not detected at IS/im, we used consecutive pairs of lO/im and 18/xm a CMa 
images to determine average offset positions between the two filters, and used these to place 
the aperture at the predicted position in the ISfim images of the Orion stars. In this way 
we measured upper limits to the flux in this band. 

The mid-IR measurements are hsted in Table 2. All stars were detected at lO/zm but 
only one measurement was secured at 18/im. We list 3-a upper limits in Table 2 for the non 
detections. 



3. Results 

3.1. Environment of the sample 

The galactic coordinates of stars in the photometric sample are plotted in Figure 1 com- 
pared to the entire sample in Paper I. The stars arc projected against the map of integrated 
^^CO emissivity from Bally et al. (1987). It can be seen that our sample is representa- 
tive of that in Paper I. It includes the CTTS and a number of WTTS in the "clump" at 
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{l = 202°,6 = -18°) in Ori OBla. It also includes a number of CTTS in Ori OBlb, located 
mostly in the inside of the molecular ring of higher extinction (cf. Paper I). Some of these 
stars are probably associated with NGC 2023 and NGC 2024. Here we adopt the boundaries 
of Ori OBlb from Paper 1, shown in Figure 1. However, there is some uncertainty in the 
actual membership of stars located near the assumed boundary between Ori OBla and lb 
at / ~ 203.7°, b ~ —18° to —19°. In addition, with the adopted boundaries, 2 stars of 
the sample fall in Ori OBlc, in the region of high molecular gas density; however, they are 
so close to the uncertain boundary that we include them in Ori OBlb for the rest of the 
analysis. 

3.2. Stelleir Properties 

We took the spectral types of the sample stars from Paper I and recalculate here stellar 
properties from the simultaneous multiwavelength photometry to minimize effects of vari- 
abihty. We used the average distances of 330 pc for Ori OBla and 440 for Ori OBlb, as 
discussed in Paper I. The extinctions Ay in Table 3 are calculated from the V — Ic color, 
using the Cardelli, Clayton, & Mathis (1989) extinction law with Ry — 3.1. The mean 
extinctions for la, lb, Ic, and the entire sample are 0.5, 0.6, 2, and 0.6, respectively. Using 
individual extinctions, we recalculated stellar luminosities L and radii R. Location in the 
HR diagram and comparison with the Siess et al. (2000) tracks yielded stellar masses M and 
ages. These quantities, shown in Table 3, in general agree within 50% with those of Paper I. 

3.3. Color-color diagrams 

Figures 2 and 3 show photometry from Table 1 plotted in the U — V vs. V — Ic and 
U — B Ys. V — Ic diagrams, as well as the location of the main sequence colors taken from 
Kenyon & Hartmann (1995, KH95). CTTS and WTTS are shown separately, and the CTTS 
in Ori OBla are marked. Note that the two observations for stars CVSO 35, 146, and 190 
are plotted as different points. 

The WTTS in the sample are consistent with main sequence colors, reddened by ~ 0.5- 
0.6 magnitudes of extinction, as expected from the average Ay of the sample. The CTTS 
have similar ^ — Jc as the WTTS, but show large excesses relative to the WTTS in U — V 
and U — B. These excesses are expected to arise from emission of the accretion shock formed 
as material accreting from the disk strikes the stellar surface (Hartigan et al. 1991; Hartigan, 
Edwards, & Ghandour 1995; Calvet & GuUbring 1998). 
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In Figures 2 and 3 we also plot a sample of stars in Taurus taken from KH95, equally 
separated in WTTS and CTTS. The age of the Taurus association has been analyzed in 
detail by Hartmann (2003), who showed that most of the stars have ages of ~ 1 - 2 Myr. 
With the ages for the subassociations found in Paper I, i.e. ~ 3 - 5 Myr for Ori OBlb and ~ 
7-10 Myr for OBla, the Taurus population is younger than Ori OBlb by ~ 2 - 3 Myr, and 
than Ori OBla by ~ 6 - 8 Myr, so it provides an initial reference point to study evolutionary 
effects. 

The Orion OBI and Taurus stars encompass a similar range of V — Ic in Figures 2 and 
3, indicating that we are probing the same range of stellar masses. The samples of Ori OBlb 
and Taurus show a similar range in U — V and U — B excesses relative to photospheric 
colors; the lowest excesses correspond to stars in Ori OBla, some of them compatible with 
measurements for the WTTS. 

In Figure 4 we show the Ori OBI sample in the U — B vs. H — K diagram, indicating the 
standard main sequence colors with values from KH95. We have used the transformations 
from Carpenter (2001) to convert the 2MASS H — Kg colors to the CIT system. A sample 
of Taurus stars, with colors from KH95 is also shown. It is apparent that the excesses in 
H — K relative to photospheric colors in Ori OBlb are smaller than those found in the Taurus 
sample. Moreover, the excesses in if — in Ori OBla, the oldest of the three populations 
shown, are even smaller than in Ori OBlb. 

This effect is more apparent in the U — B vs. K — L diagram, shown in Figure 5, using 
the X - L colors of the CTTS in Ori OBI in Table 1. We use here Ks for the Orion OBI 
stars. Colors for the Taurus stars and for the main sequence are from KH95. The color 
excesses in K — L relative to the main sequence in the Ori OBI stars are clearly smaller than 
those in Taurus. 

Finally, we locate the 6 stars observed with GEMINI/OSCIR in the iC - vs. K - L 
diagram in Figure 6, and compare it to the Taurus sample. The Ori OBI stars are located 
at the low end of the Taurus CTTS distribution, and even tend to populate the clear gap 
between CTTS and TTS noticed by KH95 (the three Taurus stars inside the gap are binaries, 
KH95). 

These comparisons clearly indicate a decrease of IR emission at < lO/xm with age. A 
significant decrease of H — K and K ~ L with with age is also found in a sample of clusters 
and groupings by Hillenbrand, Carpenter, & Meyer (2004, in preparation); Sicilia-Aguilar 
et al. (2004, in preparation) observe a similar near-IR flux deficit relative to Taurus in 
Spitzer/IRAC colors of the ~ 5 Myr cluster Tr 37. IR emission at < 10/xm arises in the 
inner disks in T Tauri stars, < few AU, which are sufficiently hot to produce such emission 
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(Meyer et al. 1997; D'Alessio et al. 1998; 1999; 2001). Thus, the observed decrease of fluxes 
with age must be related to physical phenomena occurring in these inner regions of the disks 
as a result of disk evolution. We will explore possible implications of these results in §4. 

3.4. Spectral energy distributions 

Analyses of color-color diagrams in §3.3 indicate a clear difference between emission of 
the disks in Ori OBI and in Taurus; disks in the Orion OBI stars emit less flux at near- 
infrared wavelengths than typical disks in Taurus. The comparison is seen clearly in Figure 
7, where we plot the spectral energy distributions (SEDs) of the stars observed with Gemini 
OSCIR, including the measurement and upper limits at 18 fim. from Table 2. Optical and 
near-IR fluxes have been calculated from magnitudes in Table 1 and corrected for reddening 
in Table 3. These SEDs are compared with the median SED of CTTS in Taurus (D'Alessio 
et al. 2001), scaled to the flux at H of each star. The error bars in the median show the 
flrst quartiles of the distribution (i.e., 50% of the Taurus stars fall within the error bars). 

The SEDs of the Orion stars show excesses relative to the photosphere, consistent with 
the fact that the sample observed with Gemini was selected on the basis of being bright at 
K and have large H — K excesses. However, even with this selection bias, the Ori OBI disks 
are fainter than ~ 50 - 75 % of the Taurus sample in the mid-IR. 

3.5. Mciss accretion rates 

We interpret the excesses in U — B and U — V relative to photospheric fluxes in the 
CTTS in Ori OBla and lb as due to emission from the accretion shock on the stellar surface 
(Hartigan et al. 1991; Hartigan, Edwards, & Ghandour 1995; GuUbring et al. 1998; Calvet 
& GuUbring 1998). In the framework of this interpretation, the excess luminosity above 
photospheric colors is a measure of the accretion luminosity released as matter falls from the 
disk onto the star, following the magnetic field lines that disrupt the disk (cf. Hartmann 
1988). We used the de- reddened U photometry and spectral types to obtain the excess 
luminosity in U above the intrinsic photospheric fluxes, Lu; this luminosity was used to 
estimate the accretion luminosity Lace, using the calibration of GuUbring et al. (1998). 
These estimates are listed in Table 3. The largest uncertainty in the determination of Lu, 
and thus Lace, comes from uncertainties in the determination of the extinction Ay, and it 
amounts to ~ 30 — 40%. 

Using 

Lace ^ 0.8GMM / -R*, (assuming a radius of ~ 5-R* for the magnetospheric radius; 
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Calvet & GuUbring 1998), with the mass and radius determined from the location in the 

HR diagram, we can infer the mass accretion rate M of the objects. These rates are hsted 
in Table 3. In addition to the uncertainties in Lacc^ the formal uncertainty in M includes 
errors in radius and mass, which arc of order 15-20%; thus, the mass accretion rates are 
constrained within a factor of 2. Given the uncertainties in the original calibration between 
Lu and Lacci this uncertainty can be a factor of 2 larger (GuUbring et al. 1998). As the 
errors are uncorrelated, the total error amounts to a factor of 3. 

In Figure 8 wc show histograms with the number of objects per bin of log M, each of 
size 0.5, for the three populations, Ori OBla at the top, lb at the middle, and Taurus at 
the bottom. The mass accretion rates in Taurus are taken from GuUbring ct al. (1998) and 
Hartmann et al. (1998). These measurements of M were obtained from Lace using masses 
estimated with the D'Antona & MaziteUi (1994) evolutionary tracks; we have recalculated 
them using masses for the Taurus stars estimated with the Siess et al. (2000) evolutionary 
tracks, for consistency with the work presented here. The actual values of M in bins with 
M < 10^^^MQyr~^ are highly uncertain. At these levels, emission due to stellar surface 
activity becomes an important or dominant contributor to the U excess, as comparison with 
the WTTS in the sample shows. In these cases, the determined values are upper limits to 
the true Ms. 

The age of the three populations shown increases from bottom to top. In all three 
cases, there is a large spread in the ranges of values of M. However, as time proceeds, the 
number of rapid accretors decreases, only the low accretors remain in la. One "continuum" 
star is found in Ori OB lb (see Table 3); we did not have enough spectral resolution to 
distinguish spectral features and thus determine its spectral type and derived properties. 
This continuum star does not appear in Figure 8; if its properties were similar to those of 
continuum stars in Taurus, then it could have M ~ few lO^^M© yr~^ (GuUbring et al. 1998; 
Calvet & GuUbring 1998; GuUbring et al. 2000). However, even including this star, the 
relative number of rapid accretors in Ori OB lb would be smaller than in Taurus. 

We compared the distributions of log M in the three regions using a two-sided KS test. 
The results indicated that the probability that the Ori OBlb and Taurus samples are drawn 
from the same underlying distribution is only 3%; the same is true for Ori OBlb vs. Ori 
OBla. For a further test, we have added the continuum star in Ori OBlb to the distribution, 
assuming a a mass accretion rate of lO^^M© yr~^; even in this case, the probability that the 
Ori OBlb sample is drawn from the same distribution as Taurus is only marginally larger, 
9%. Thus we conclude that the distributions of accretion rates are significantly different 
between the three regions. 

Figure 9 plot the derived mass accretion rates and ages on CTTS in Ori OBla and lb 
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in Table 3. We also plot in Figure 9 similar data for Taurus, p Ophiucus and Chameleon 
(Hartmann et al. 1998), and the TW Hya association (TWA, MuzeroUe et al. 2000, 2002). 
All masses and mass accretion rates of these latter objects have been re-calibrated using the 
Siess et al. (2000) evolutionary tracks. The mass accretion rates of the disks in Ori OBI 
are consistent with those in other associations, and show a clear decrease with age. The star 
CVSO 41 in Ori OB la has an apparent age of 88 Myr (Table 3) and still shows a fairly high 
value of M. The spectral type of this star is K2, so it could have a large age uncertainty. 
It is known that pre-main-sequence stars in this spectral type range tend to have higher 
apparent ages than cooler objects (Hillenbrand 1997); this behavior has been attributed to 
uncertainties in the birthline location (Hartmann 2003). 

The solid hne in Figure 9 shows the decay of M with time expected from viscous evolu- 
tion for a disk with initial mass O.2M0 and viscous parameter a — 0.01, from Hartmann et 
al. (1998). The new values of M help confirm a decrease of M with age qualitatively consis- 
tent with viscous disk evolution. However, viscous evolution cannot explain the decreasing 
fraction of accreting objects with age (Muzerolle et al. 2000; Briceho et al. 2004); other 
factors (i.e., inner disk clearing associated with planet formation, Calvet et al. 2002; Rice et 
al. 2003; D'Alessio et al. 2004) must play a role in slowing accretion onto the central star. 

4. Discussion 

In our comparison of disk properties in three populations with ages ~ 1-2 Myr (Taurus), 
3-5 Myr (Ori OBlb), and 7-10 Myr (Ori OBla), we have found a clear decrease of disk 
emission with age. This decrease is apparent at all wavelengths with expected substantial 
disk contribution to the emission, from H to 18/xm. 

The disk emission in the IR depends on the temperature of the dust. The two main 
heating mechanisms active on disks around young stars are viscous dissipation and stellar 
irradiation (Kenyon & Hartmann 1987, KH87; D'Alessio et al. 1998, 1999, 2001, DOl). 
Viscous dissipation depends on the mass accretion rate throughout the disk; for the range of 
stellar parameters of stars in Ori OBI, and for M < IO^^MqUT^^, which characterize these 
stars (§3.5), viscous heating is important only near the midplane of the innermost disk, < 1 
AU (D'Alessio et al. 1999, 2001). Thus, absorption of stellar radiation by dust in the disks 
around the Ori OBI stars mostly determines their temperature and emission. 

One important development in recent years is the recognition that the near-IR fluxes 
in accretion disks around pre-main sequence stars are dominated by emission from a sharp 
transition between dust and gas at the dust destruction radius. At this transition, the 



- 12 - 



optically thick dust disk is high enough to form a "wall" , which is frontally illuminated by 
the star. Opacities in the inner gas disk are due to molecules and are much lower than 
the dust opacities for typically low mas accretion rates; thus, the gas disk is optically thin, 
enabling frontal illumination of the wall (Muzcrollc et al. 2004). The existence of this 
wall was proposed for Herbig Ae/Be stars (Natta ct al. 2001; Tuthill, Monnier & Danchi 
2001; DuUemond, Dominik, & Natta et al. 2001) to explain the peculiar "bump" at ~ S/xm 
observed in the SEDs of these stars (Hillenbrand et al. 1992). Interferometric measurements 
of inner disk structures in Herbig Ae/Be stars confirmed this prediction (Millan-Gabet et al. 
1999; Monnier & Millan-Gabet 2002). A similar wall emission was then reported in the low 
mass CTTS by MuzeroUe et al. (2003). In contrast to their higher mass counterparts, CTTS 
can have accretion luminosities of the same order as the stellar luminosities, thus increasing 
the total energy output and moving outwards the disk location where the temperature is 
equal to the dust destruction temperature for objects with high values of Lace- The rapid 
accretors are therefore expected to have the wall located at larger radius than objects with low 
M, increasing the total emitting area of the wall and thus the near-IR excess (MuzeroUe et al. 
2003; D'Alessio et al. 2003). This prediction has been confirmed by the Keck interferometer 
measurement of an inner structure at the expected location of the dust destruction radius 
in a high accretion CTTS (Colavita et al. 2003; D'Alessio et al. 2003). 

According to this model, the near-IR excess depends on the strength of the wall emission 
relative to photospheric fluxes. This emission is approximately that of a black-body at the 
dust destruction temperature, ~ 1400K (MuzeroUe et al. 2003), times the sohd angle 
subtended by the wall. This sohd angle is essentially that of a cyhnder of given height, and 
radius equal to the dust destruction radius, with an inclination dependence which makes the 
emission maximum at intermediate lines of sight (60° - 80°) where the wall has maximum 
exposure (DuUemond et al. 2001; D'Alessio et al. 2003). The dust destruction radius 
is essentially oc (L^ + Lacc)^^"^ /Tj (D'Alessio ct al. 2003). The height of the wall 
is estimated as the height where the radial optical depth from the star becomes ~ 1 (see 
MuzeroUe et al. 2004). 

Wc can examine the decrease of near-IR excesses with age in terms of these ideas. As 
they age, pre-main sequence low mass stars move down along Hayashi tracks essentially at 
constant T^ff while decreasing in luminosity. As M decreases. Lace becomes much smaller 
than L*, so the dust destruction radius becomes fixed relative to the stellar radius, Ra/ R* oc 
L\^'^ /T^R^ ~ {^/TdY ~ const. If the height of the wall is a fixed multiple of the disk scale 
height at Ra, which in turn is typically ~ O.lii!^, then the emitting area of the wall, oc RdZ^, 
and thus its total energy output, would decrease in the same proportion as the stellar output, 
namely oc R^. In this case, no decrease of near-IR excess relative to the photosphere with 
time would be expected. Since Rd/R* ~ const, the most natural way to explain the observed 
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decrease is that the height of the wall decreases with time. For this to happen, the opacity 
at the upper levels of the wall has to decrease, which implies that dust grains at the top of 
the wall significantly grow and/or settle towards the midplane as age increases. 

Grain growth and/or settling towards the midplane in the disk regions outside the wall 
also result in low mid-IR fluxes. The disk is heated by stellar + accretion energy captured 
at the disk surface. Thus, heating of the disk depends on the shape of this surface; the more 
"flared" the surface is, the more energy it can capture (Kenyon & Hartmann 1987, KH87); 
disks with highly flared surfaces are hotter and thus emit more. To be more speciflc, the 
height of the surface Zs at a given radius is given by the condition that the radial optical Since 
the density decreases with vertical height above the midplane, Zs decreases as the opacity 
of the dust in the upper layers decreases. Small grains have high opacity at wavelengths 
characteristic of the stellar radiation, in the optical and near-IR. They produce a highly 
flared surface, which results in high energy capture, and hot and bright disks. In contrast, 
large grains have low opacity at short wavelengths. They produce flatter surfaces, which 
result in less energy capture and cooler and fainter disks (D'Alessio et al. 2001, DOl). 

Alternatively, the opacity at upper layers may be reduced because grains are settling to 
the midplane. This also results in flatter surfaces and fainter disks in the infrared (DuUemond 

& Dominik 2004; D'Alessio et al. 2004, in preparation). Both growth and settling are 
expected to occur together. Dust particles, affected by gas drag, collide and stick together; 
as the particles grow, they settle towards the midplane, falling under the effect of gravity 
(Wcidcnschilling & Cuzzi 1993; Wcidcnschilling 1997). These theories have many unknown 
parameters, and so far predict time scales for complete settling < 1 Myr; however, the 
timescales depend upon the magnitude of turbulence in disks, which could retard or prevent 
setthng. 

A number of claims of grain growth in disks have been made over the years, starting 
with the interpretation of slopes of the fluxes in the millimeter, which are expected to become 
flatter as the grain size increases (Beckwith & Sargent 1991). Assuming a dust mixture and 
a dust size distribution of the form n{a) oc a'^ between minimum and maximum sizes amin 
and Umax, DOl flnd that a much better flt to the median SED of Taurus can be achieved with 
dust particle growth to a maximum size amax of 1 n:m:i (p=3.5) than with interstellar grains, 
with ttmax — 0.25/im. Grain growth and settling has been claimed to explain the SEDs of 
some CTTS in Taurus. One well studied case is that of the CTTS LkCa 15, which has a 
SED much flatter than the Taurus median in the IR and high emission at submillimeter and 
millimeter wavelengths (Qi et al. 2003); the interpretation of this SED requires signiflcant 
grain growth and settling (Chiang et al. 2001; D'Alessio et al. 2003, 2004; Bergin et al. 
2004). In Figures 2, 3, 4, 5, and 6 we have marked the location of LkCa 15; disks in Ori 
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OBl have similar colors to those of LkCa 15 in all bands, suggesting that significant dust 
evolution has also taken place in these older disks. In Figure 10 we compare the SED of 
CVSO 109, the only object for which we have detections up to ISfim, to fluxes for LkCa 
15 from KH95 and to the Taurus median, both scaled at the H of CVSO 109. Note that 
LkCa 15 is a K5 star (KH95), so it is brighter in the optical than CVSO 109, which is an 
MO (Table 3). The SEDs of LkCa 15 and CVSO 109 are significantly lower than the Taurus 
median, although the flux deficit is much larger in CVSO 109. 

To try to quantify the degree of dust evolution in CVSO 109, we have calculated models 
with uniformly distributed gas and dust following the procedures of DOl. We have used 
stellar and accretion parameters from Table 3, and have varied the parameters ttmax and p 
which characterize the dust size distribution. Figure 10 shows our best model. To fit the 
observed fluxes, we require Umax = 10 cm; this is higher than the maximum sizes required 
to explain the median SED in Taurus, Umax ~ 1 mm (DOl), or the SED of the ~ 10 Myr old 
star TW Hya, a^ax ~ 1 cm (Calvet et al. 2002; Wilner et al. 2004). Moreover, we require 
p = 2.5, a value consistent with significant dust coagulation (Miyake & Nakagawa 1995); 
modeling of CTTS disks can generally be carried out successfully with the assumption of 
p = 3.5, a value characteristic of the interstellar medium (cf. DOl, Chiang et al. 2001). 

Our modeling is only indicative, in that we have assumed that the dust is uniformly 
distributed with the gas. Even with this assumption, our comparison implies that the disk of 
CVSO 109 is essentially a flat disk, with little if any flaring. Thus, our modeling indicates that 
the dust in the disk of CVSO 109 has strongly settled towards the midplane, as expected from 
theories of dust evolution in protoplanetary disks (Weidenschilling 1997). Models including 
settling in addition to observations at longer wavelengths are required to actually constrain 
the size distribution of the solids in the midplane. 

A possible alternative to consider is that disks in Ori OBI are smaller than their Taurus 
counterparts because their outer parts have been photoevaporated by the external ionizing 
fields of the high mass stars in the associations. A smaller disk would have less emission 
at long wavelengths since the low temperature regions would be removed. Models for the 
viscous evolution of disks subject to the effects of photoevaporating radiation fields from 
the accretion flows and external radiation fields have been calculated by Clarke, Gendrin, & 
Sotomayor (2001) and Matsuyama, Johnston & Hartmann (2003). Since the distance from 
the ionizing stars are of the order of a few pc in OB associations, FUV external photons 
are expected to dominate over EUV photons in the photoevaporation of disks (Johnstone et 
al. 1998). For this case, Matsuyama et al. (2003) find that that the outer portions of the 
disk could be removed on time scales which range from 1 to several Myr, depending on the 
viscosity parameter. After significant photoevaporation, the disk edge is fixed by the FUV 
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gravitational radius, which for parameters as those of the Ori OBlb stellar sample is ~ 30 
AU. To determine the effects of a reduction of disk radius on the SED, and in particular to 
see if a smaller disk radius can explain the low mid-infrared fluxes without extreme grain 
growth, we have calculated disk models with ttmax = 1 and p = 3.5, characteristic of the 
Taurus population, for disk radii 100 AU, 50 AU, and 10 AU. These SEDs are shown in Figure 
11. Models are calculated at a high inclination to the line of sight (72°) to maximize the 
contrast between the star and the disk. Even small disks produce too much long wavelength 
flux compared to the the observations, if a^ax — 1 mm ^ in Taurus. This implies again that 
the disk has to be essentially flat, which require that grains grow to larger sizes and/or have 
settled. 



5. Summary and conclusions 

We have detected clear signatures of disk evolution in samples of accreting stars in the 
Ori OBla and lb associations, with ages 7-10 Myr and 3-5 Myr, respectively. The mass 
accretion rates of these stars are consistent with viscous evolution of accretion disks. In 
addition, we find a significant overall decrease of infrared emission with age. Comparison 
with disk models indicates that grain growth and settling, processes expected from dust 
evolution, are responsible for this decrease. Future observations of larger samples of Ori 
OBI stars with SPITZER will provide more detail on the nature of the dust evolution 
process. 
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Table 1. OPTICAL AND NEAR INFRARED PHOTOMETRY 
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Table 1 — Continued 
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157 


17.87 


0.70 


-0.98 


1.14 


2.09 


14.01 


0.59 


0.35 




164 


16.71 


1.55 


0.30 


1.31 


3.03 


11.80 


0.69 


0.32 




165 


13.68 


1.12 


0.45 


0.76 


1.48 


11.07 


0.83 


0.41 


9.24 


168 


17.33 


1.24 


-0.25 


1.27 


2.82 


12.78 


0.60 


0.28 




171 


14.84 


1.46 


0.64 


0.84 


1.69 


11.96 


0.70 


0.30 




176 


15.65 


1.21 


-0.35 


1.14 


2.46 


11.71 


0.91 


0.50 


9.84 


177 


15.62 


0.71 




0.96 


1.92 


12.30 


0.77 


0.33 


10.98 
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Table 1 — Continued 



cvso 


V 


B-V 


U-V 


V-R 


V-I 


J 


J-H 


H-K, 


L 


178 


15.66 


1.48 


0.64 


1.11 


2.28 


11.47 


0.77 


0.34 




181 


15.63 


1.22 




0.87 


2.00 


12.30 


0.74 


0.19 




184 


15.34 


1.86 




1.41 


2.75 


10.40 


1.27 


0.91 




185 


15.28 


1.61 




1.08 


2.18 


10.88 


0.88 


0.49 


8.77 


190 


14.04 


1.10 


0.28 


0.83 


1.72 


11.04 


1.08 


0.77 


8.49 




14.76 








1.89 


11.04 


1.08 


0.77 


8.49 


192 


15.99 


1.26 


-0.14 


1.17 


2.37 


11.56 


1.01 


0.55 


99.00 


193 


15.15 


1.29 


0.04 


1.15 


2.29 


11.34 


1.22 


0.66 


99.00 



Table 2. OSCIR PHOTOMETRY 



CVSO 


N 


errN 


[18/xm] 


err[18/xm] 


log \Fx{N) 


log AFA(18/xm) 












log (erg sec~^ cm~^) 


log (erg sec"^ cm~^) 


47 


9.418 


0.081 


> 7.627 


0.383 


-11.779 ± 0.03 


< -11.654 


58 


8.448 


0.388 


> 7.517 


0.444 


-11.391± 0.15 


< -11.576 


90 


8.294 


0.282 


> 9.533 


2.773 


-11.330± 0.12 


< -11.792 


107 


8.109 


0.276 


> 9.092 


1.866 


-11.256± 0.12 


< -11.764 


109 


7.935 


0.120 


6.437 


0.165 


-11.186± 0.05 


-11.492± 0.1 


165 


6.921 


0.093 


> 6.983 


0.279 


-10.781± 0.05 


< -11.710 
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Table 3. Stellar and Accretion Properties 



cvso 


SpT 




L 


R 


M 


Av 


Age 


Lace 


M 


Assoc 


Type 






K 


Lq 


Rq 


Mq 




Myr 


Lq 


lO-8M0j/r--i 






1 


MS 


3470 


0.18 


1.17 


0.33 


0.81 


4.02 


0.01 


0.09 


la 


CTTS 


3 


M3 


3470 


0.16 


1.11 


0.S3 


0.00 


4.68 


0.00 


0.00 


la 


WTTS 


5 


M2 


3580 


0.21 


1.20 


0.38 


0.22 


4.12 


0.00 


0.01 


la 


WTTS 


10 


MO 


3850 


0.27 


1.17 


0.58 


0.27 


7.20 


0.00 


0.00 


la 


WTTS 


11 


M2 


3580 


0.33 


1.49 


0.39 


0.46 


2.51 


0.00 


0.02 


la 


WTTS 


15 


Ml 


3720 


0.21 


1.11 


0.47 


0.39 


6.42 


0.00 


0.03 


la 


WTTS 


17 


MS 


3470 


0.19 


1.20 


0.S3 


0.00 


3.63 


0.00 


0.06 


la 


WTTS 


19 


MS 


3470 


0.34 


1.60 


0.34 


0.10 


2.23 






la 


WTTS 


22 


M2 


3580 


0.18 


1.12 


0.38 


0.36 


5.18 


0.00 


0.01 


la 


WTTS 


24 


M2 


3580 


0.17 


1.06 


0.38 


0.27 


5.97 


0.00 


0.01 


la 


WTTS 


25 


MO 


3850 


0.24 


1.11 


0.58 


0.27 


8.09 






la 


WTTS 


29 


M2 


3580 


0.29 


1.40 


0.39 


0.94 


2.73 






la 


WTTS 


34 


MO 


3850 


0.29 


1.21 


0.58 


0.41 


6.61 


0.00 


0.01 


la 


WTTS 


34 


MO 


3850 


0.19 


0.98 


0.59 


0.41 


10.55 


0.00 


0.01 


la 


WTTS 


35 


K7 


4060 


0.72 


1.72 


0.77 


0.85 


2.88 


0.02 


0.17 


la 


CTTS 




K7 


4060 


0.71 


1.70 


0.77 


0.80 


2.90 


0.02 


0.13 


la 


CTTS 


36 


M3 


3470 


0.21 


1.28 


0.34 


0.00 


2.95 


0.03 


0.36 


la 


WTTS 


38 


M2 


3580 


0.17 


1.06 


0.38 


0.31 


5.88 


0.00 


0.00 


la 


WTTS 


38 


M2 


3580 


0.19 


1.13 


0.38 


0.77 


4.98 


0.00 


0.02 


la 


WTTS 


39 


M2 


3580 


0.30 


1.41 


0.39 


0.74 


2.70 


0.00 


0.04 


la 


WTTS 


39 


M2 


3580 


0.30 


1.42 


0.39 


0.77 


2.68 


0.01 


0.08 


la 


WTTS 


40 


MO 


3850 


0.24 


1.10 


0.58 


0.02 


8.18 


0.00 


0.00 


la 


CTTS 


41 


K2 


4900 


0.37 


0.84 


0.89 


1.30 


88.18 


0.06 


0.19 


la 


CTTS 


42 


M2 


3580 


0.17 


1.08 


0.38 


0.02 


5.71 


0.00 


0.00 


la 


WTTS 


43 


M2 


3580 


0.25 


1.29 


0.39 


0.77 


3.03 


0.00 


0.02 


la 


WTTS 


44 


K6 


4205 


0.27 


0.98 


0.82 


0.00 


22.97 


0.00 


0.01 


la 


WTTS 


46 


MO 


3850 


0.36 


1.S4 


0.58 


0.32 


4.74 


0.00 


0.03 


la 


WTTS 


46 


MO 


3850 


0.35 


1.S4 


0.58 


0.29 


4.80 


0.00 


0.02 


la 


WTTS 


47 


K5 


4350 


0.97 


1.73 


1.10 


1.21 


4.96 


0.01 


0.04 


la 


CTTS 


48 


K4 


4590 


1.25 


1.77 


1.28 


0.71 


6.60 


0.05 


0.22 


lb 


CTTS 


54 


MO 


3850 


0.17 


0.92 


0.58 


0.36 


14.57 


0.00 


0.01 


la 


WTTS 


55 


M2 


3580 


0.37 


1.58 


0.39 


0.31 


2.30 


0.00 


0.01 


la 


WTTS 


55 


M2 


3580 


0.37 


1.58 


0.39 


0.29 


2.31 


0.00 


0.01 


la 


WTTS 


57 


MO 


3850 


0.20 


1.00 


0.59 


0.27 


9.95 


0.00 


0.00 


la 


CTTS 


58 


K7 


4060 


0.61 


1.58 


0.78 


0.12 


4.02 


0.07 


0.45 


lb 


CTTS 


72 


M6 


3050 


0.04 


0.74 


0.11 


0.00 


6.20 


0.00 


0.01 


lb 


CTTS 


75 


Ml 


3720 


0.23 


1.16 


0.47 


0.00 


5.82 


0.00 


0.04 


lb 


CTTS 


77 


K6 


4205 


0.28 


0.99 


0.82 


1.88 


21.91 






la 


WTTS 


81 


MS 


3470 


0.45 


1.86 


0.34 


0.17 


1.76 


0.00 


0.02 


lb 


WTTS 


90 


K7 


4060 


0.55 


1.50 


0.79 


0.00 


5.04 


0.29 


1.77 


lb 


CTTS 


104 


K7 


4060 


0.76 


1.76 


0.77 


0.00 


2.76 


0.10 


0.75 


lb 


CTTS 


107 


K7 


4060 


1.02 


2.04 


0.77 


0.32 


2.12 


0.13 


1.09 


lb 


CTTS 


109 


MO 


3850 


1.48 


2.73 


0.56 


0.00 


0.80 


0.16 


2.52 


lb 


CTTS 


121 


K3 


4730 


1.31 


1.70 


1.32 


1.49 


8.06 


0.40 


1.67 


lb 


CTTS 
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Table 3 — Continued 



cvso 


SpT 




L 


R 


M 


Av 


Age 


Lace 


M 


Assoc 


Type 






K 


Lp, 




M';.-: 




Myr 


L,:;, 








124 


M3 


3470 


0.51 


1.97 


0.34 


0.00 


1.58 


0.00 


0.02 


lb 


WTTS 


126 


MO 


3850 


0.16 


0.89 


0.58 


0.12 


16.16 


0.00 


0.00 


la 


WTTS 


133 


K6 


4205 


0.32 


1.06 


0.84 


0.27 


18.14 


0.00 


0.01 


la 


WTTS 


143 


Ml 


3720 


0.82 


2.18 


0.47 


0.55 


1.16 


0.02 


0.37 


lb 


CTTS 


146 


K6 


4205 


1.32 


2.17 


0.91 


0.37 


2.01 


0.11 


0.81 


lb 


CTTS 




K6 


4205 


1.23 


2.09 


0.92 


0.10 


2.19 


0.12 


0.88 


lb 


CTTS 


152 


K3 


4730 


1.21 


1.64 


1.30 


0.78 


8.65 


0.03 


0.14 


lb 


CTTS 


153 


M2 


3580 


0.23 


1.26 


0.39 


0.36 


3.43 


0.02 


0.16 


la 


WTTS 


155 


K6 


4205 


1.30 


2.14 


0.91 


0.90 


2.06 


0.07 


0.53 


lb 


CTTS 


157 


C 


















lb 


CTTS 


164 


M3 


3470 


1.06 


2.85 


0.35 


1.82 


0.45 


0.05 


1.19 


lb 


CTTS 


165 


K6 


4205 


1.48 


2.29 


0.91 


0.00 


1.72 


0.05 


0.37 


lb 


CTTS 


168 


M4 


3370 


0.25 


1.46 


0.29 


0.00 


2.49 


0.00 


0.05 


lb 


CTTS 


171 


K5 


4350 


0.82 


1.60 


1.08 


0.81 


6.41 


0.02 


0.10 


lb 


CTTS 


176 


M3 


3470 


0.71 


2.32 


0.35 


0.00 


1.06 


0.02 


0.43 


lb 


CTTS 


177 


K6 


4205 


0.64 


1.50 


0.94 


1.07 


6.21 


0.01 


0.04 


lb 


CTTS 


178 


Ml 


3720 


1.17 


2.60 


0.47 


0.77 


0.81 


0.01 


0.25 


lb 


CTTS 


181 


MO 


3850 


0.51 


1.61 


0.57 


0.48 


2.63 


0.02 


0.14 


lb 


CTTS 


184 


Kl 


5080 


9.67 


4.01 


2.35 


4.35 


2.03 






lb 


CTTS 


185 


K7 


4060 


2.54 


3.22 


0.81 


1.40 


0.79 






lb 


CTTS 


190 


K6 


4205 


1.78 


2.52 


0.90 


0.58 


1.26 


0.13 


1.17 


lb 


CTTS 




K6 


4205 


1.99 


2.66 


0.90 


0.99 


1.00 


0.14 


1.38 


lb 


CTTS 


192 


K6 


4205 


1.68 


2.44 


0.90 


2.14 


1.40 


0.36 


3.08 


Ic 


CTTS 


193 


K6 


4205 


1.96 


2.63 


0.90 


1.95 


1.02 


0.48 


4.49 


Ic 


CTTS 
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^B. ;P3^- :;^<J 



i 

Fig. 1. — Galactic coordinates of the stars in the photometric sample (big open circles) 
compared to the Orion Variability Survey sample of Paper I (CTTS, filled circles, WTTS 
filled triangles). The integrated ^"^CO emissivity from Bally et al. (1987) is shown in a 
halftone grey scale, covering the range from 0.5 to 20 Kkms~^. The isocontour for Ay — 1 
from Schlegel et al. (1998) is also shown. The dotted lines show the boundaries of Ori OBlb 
adopted in Paper I and this work. The dashed lines show the hmits of the survey in Paper 
I. 
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12 3 4 

V-Ie 

Fig. 2. — Location of the stars in the U — V vs. V — Ic diagram. We show the CTTS (filled 
large circles), WTTS (open large circles) for the Ori OBI sample, with the stars in Ori OB la 
encircled by the open rectangle. The Taurus CTTS (small filled circles) and WTTS (small 
open circles), and the sequence for dwarf standards from Kenyon & Hartmann (1995) are 
also shown, as well as the reddening vector for Ay = 1 calculated with the Cardelli et al. 
(1989) reddening law for Ry = 3.1. The star symbol corresponds to the Taurus CTTS LkCa 
15. 
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Fig. 3. — Location of the stars in the U — B vs. V — Ic diagram. Symbols as in Figure 2. 
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Fig. 4. — Location of the stars in the U — B vs. H — K diagram. Symbols as in FigTirc 2. 
The H — K colors for the Ori OBI stars are from the 2MASS survey, corrected to the CIT 
system. 
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Fig. 5. — Location of the stars in the U — B vs. K — L diagram. Symbols as in Figure 2 
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Fig. 6. — Location of the stars in the K — L vs. K — N diagram. Symbols as in Figure 2. 
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Fig. 7. — Spectral energy distributions of stars observed with GEMINI/OSCIR (filled circles) 
compared to the Taurus median (D'Alessio et al. 2001, dashed lines), normalized to the H 
magnitude of each star. The error bars in the Taurus median correspond to the first quartiles. 
The inverted triangles indicate upper limits. The photospheric fluxes in the near-IR (dotted 
lines) are constructed from colors of dwarfs from Kenyon & Hartmann (1995), normalized 
to the H magnitude of each star. 
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log dM/dt 

Fig. 8. — Distribution of mass accretion rate for Oril OBla (upper panel), Ori OBlb (middle 
panel) and Taurus (lower panel). Data from Taurus from GuUbring et al. (1998) and 
Hartmann et al. (1998), recalculated with masses estimated from the Siess et al. (2000) 
evolutionary tracks. 
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log age(yrs) 



Fig. 9. — Mass accretion rate vs. age for several associations: Ori OBla (open triangles), 
Ori OBlb (open rectangles), Taurus (filled circles), p Oph (open circles), Chamaeleon I 
(filled triangles), TWA (inverted open triangles), with data from Hartmann et al. (1998) 
and Muzerolle et al. (2000, 2001). Solid line: model for viscous evolution for an initial disk 
mass of 0.2Mq and a = 0.01 from Hartmann et al. (1998). A typical error bar is shown. 
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Fig. 10. — Spectral energy distribution of CVSO 109 (solid filled circles) compared to the 
median of Taurus (dashed lines), the photospheric fluxes (dotted lines), and the SED of 
LkCa 15 (x's). The solid line shows the SED of an irradiated accretion disk model with 
similar stellar and accretion parameters as CVSO 109, and with a dust size distribution 
characterized by a^aa; = 10 cm and exponent p = 2.5 (see text). 
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Fig. 11. — Effect of disk radius on SED. Disk models for a„iax = Imm and p = 3.5, and disk 
radius 100 AU (long-dashed line), 50 AU (short-dashed line), and 10 AU (dot-dashed line). 
Model for a^ax — 10 cm and p — 2.5 in sohd line and observations of CVSO 109 in filled 
circles. The photosphere is indicated by the dotted line. 



